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At reduced dimensionality, Coulomb interactions play a crucial role in determining device proper- 
ties. While such interactions within the same carbon nanotube have been shown to have unexpected 
properties, device integration and multi-nanotube devices require the consideration of inter-nanotube 
interactions. We present calculations of the characteristics of planar carbon nanotube transistors 
including interactions between semiconducting nanotubes and between semiconducting and metallic 
nanotubes. The results indicate that inter-tube interactions affect both the channel behavior and 
the contacts. For long channel devices, a separation of the order of the gate oxide thickness is 
necessary to eliminate inter-nanotube effects. Because of an exponential dependence of this length 
scale on dielectric constant, very high device densities are possible by using high-K dielectrics and 
embedded contacts. 

I. INTRODUCTION 



Much recent experimental and theoretical work has fo- 
cused on the electronic transport properties of nanowires, 
motivated by promises of novel electronic devices and the 
basic scientific challenges that they present. One of the 
key findings in nanowire devices is that the Coulomb 
interaction plays a crucial role in determining device 
properties. For example, electrostatics in carbon nan- 
otubes (NTs) leads to unexpected behavior of intra-tube 
p-n junctions 1 , of intra-tube Schottky junctions 2 , and of 
Fermi level pinning at contacts between NTs and metals 3 . 
The importance of m£ra-tube Coulomb interactions nat- 
urally leads to the question of how inter-tube interac- 
tions might influence device behavior. Answering this 
question is gaining increasing importance as experimen- 
tal devices based on multiple nanotubes are becoming 
more common 4-7 [an example is shown in Fig. 1(a)], and 
also to address the question of device packing density. 

In this paper, we present self-consistent calculations 
based on the non-equilibrium Green's function technique 
for planar carbon nanotube transistors containing multi- 
ple parallel NTs. By analyzing the variations of the tran- 
sistor characteristics with NT separation, we establish a 
length scale below which inter-tube interactions become 
important. While for small channel lengths this length 
scale depends on the channel length, for long channel de- 
vices it becomes independent of the channel length and 
is essentially determined by the gate oxide thickness. Im- 
portantly, this length scale depends exponentially on the 
dielectric constant of the medium surrounding the NTs, 
and can be substantially reduced by embedding the NTs 
in a high-K dielectric in the channel and using embedded 
contacts. 



II. COMPUTATIONAL APPROACH 

Figure 1 (b-e) shows the device under consideration: in- 
finitely long single-wall zigzag NTs laying on metal elec- 
trodes at their two ends and on a dielectric in the channel 
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FIG. 1. Panel (a) shows a SEM image of a carbon nan- 
otube device containing multiple nanotubes. Panel (b) shows 
a sketch of the device used in our calculations, with the dis- 
tance L indicating the channel length and d the NT separa- 
tion. Panel (c) shows a side view of the device, while panels 
(d) and (e) show top views of the unit cells used to study 
semiconducting-semiconducting and semiconducting-metallic 
interactions, respectively. 

the NT and the metal or dielectric, and assume that 
the NT cross-section remains circular, a good approx- 
imation for the NTs considered here 8 . The electronic 
properties of the NT are described using a tight-binding 
approach with one it orbital per carbon atom, and a cou- 
pling 7 = 2.5 eV between nearest-neighbor atoms. Semi- 
conducting tubes are (17,0) NTs (band gap of 0.55 eV, 
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diameter of 1.33 nm) while metallic tubes are (18,0) NTs 
(diameter of 1.41 nm). We take the NT midgap as the 
energy reference level, and use a metal workfunction of 
5.5 eV, putting the metal Fermi level 1 eV below the NT 
midgap. The gate insulator dielectric constant is that of 
Si0 2 . 

To apply the non-equilibrium Green's function 
formalism 9,10 to this system, we divide the NT in princi- 
ple layers, with each layer corresponding to a ring of the 
zigzag NT. The main quantity of interest is the Green's 
function G R , from which one can obtain the transmission 
probability T(E) and the zero bias conductance 



C = 



2e7 



T(E) 



dfjE) 
dE 



dE 



(1) 



where / is the Fermi function. G R is calculated by solv- 
ing the matrix equation 



G R = [{E -eU)I-H - \Z R ] 



(2) 



where Ho is the tight-binding Hamiltonian for the iso- 
lated NT and U is the electrostatic potential evaluated 
at the position of each layer. The functions T, R and 
represent the interaction of the scattering region with the 
semi-infinite NT leads. In our tight-binding represen- 
tation, the Hamiltonian matrix elements for layer I are 

rr2Z,2Z-l rr2/-l,2i o„, „„„ /" 7rJ \ tt21,21 + 1 tj21+1,21 

H =H = 27 cos (^J, H -H - 
7 where M is the number of atoms around a NT ring ( 
M = 17 or 18 here) and J = 1,...,M labels each of 
the NT bands. In this representation, we assume that 
the electrostatic potential on every atom of a ring is the 
same, a good approximation in this case since U varies 
slowly over the NT diameter. 

The electrostatic potential is calculated by solving 
Poisson's equation in three-dimensional coordinates on 
a variable grid with the source charge on the NT, and 
with boundary conditions at the source, drain, and gate 
surfaces, and at the vacuum/dielectric interface. This 
is done by combining a finite difference approach in the 
directions parallel to the NT and perpendicular to the 
substrate with a fast fourier transform in the direction 
perpendicular to the NTs. Once the three-dimensional 
electrostatic potential is obtained, the value for U on 
each ring is taken as the average of the electrostatic po- 
tentials on the points on the NT nearest and furthest 
from the substrate, at the position of each ring. 

To obtain the charge density, we note that our tight- 
binding technique provides the total charge on each 
layer of the NT, which needs to be spatially distributed. 
We assume a uniform distribution of the charge in the 
azimuthal direction, and spatially distribute the total 
charge in the radial and axial directions with a Gaussian 
smearing function. The three-dimensional charge density 
is then given by 



a(r,z,, 
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where g(z—zi, r—R) = (47r 2 i?cr z oy) 1 exp \—(z — z/) 2 /2<r 2 ] 
exp [— (r — i?) 2 /2<r 2 ] with R the tube radius, zi the po- 
sition of ring Z, and <r z and oy the smearing lengths in 
the axial and radial directions respectively (this expres- 
sion for g is valid when R a r , and we used values of 
<j z = 0.14 nm and oy = 0.06 nm). Our overall proce- 
dure is to solve the coupled set of nonlinear Eqs (2 — 3) 
self-consistently for a given gate-source voltage at zero 
drain-source bias. Figure 1(e) shows the unit cells for 
our calculations to study semiconducting-semiconducting 
and semiconducting-metallic interactions. 
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FIG. 2. Panels (a) and (b) show band diagrams for d = 25.6 
nm at gate-source voltages of -4 and 8 volts, respectively. 
Panel (c) shows the conductance of a (17,0) NT transistor 
as a function of gate-source voltage, for different NT separa- 
tions. The inset shows collapse of the data after rescaling and 
shifting of the gate voltage. 



III. RESULTS AND DISCUSSION 



We first discuss interactions between semiconducting 
NTs. The conductance of the NT transistor as a function 
of the gate-source voltage is plotted in Fig. 2 for differ- 
ent values of the separation between the NTs. Clearly, 
as the NT separation is reduced, the NTFET character- 
istics degrade; in the range of gate voltages studied, the 
current cannot be completely turned on for the smallest 
tube separation, while the behavior in the OFF state is 
essentially independent of tube separation. This behav- 
ior arises because of the charge induced on the NTs by 
the gate voltage, and by interactions between NTs at the 
contact. As shown in Fig. 2a, the OFF state corresponds 
to the Fermi level in the middle of the NT bandgap in 
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the channel, with little charge on the NTs, so Coulomb 
interactions in the channel between neighboring NTs are 
negligible. The ON state however consists of making the 
NT p-type in the channel by raising the valence band 
above the Fermi level with the gate voltage (Fig. 2b); 
this creates a large positive charge on the NT which in- 
teracts with the charge on neighboring NTs, and lowers 
the electron energy. This pushes the bands down, and 
reduces the conductance. 

While interactions in the channel play an important 
role in determining device behavior, NT-NT interactions 
also affects the contacts. Because of the mismatch be- 
tween the metal and NT workfunctions, contacts between 
NTs and metals are governed by charge transfer between 
the metal and the NT, effectively doping the NT p-type 
at the contact, with the metal Fermi level just below the 
NT valence band edge in the contact 11 . As the separa- 
tion between NTs decreases, the Coulomb interactions 
between the positive charge at the contacts leads to a 
lowering of the bands; at small enough separations, the 
valence band in the contact can be pushed below the 
metal Fermi level. Thus inter-tube interactions at the 
contacts can lead to modifications of the NT/metal con- 
tacts from ohmic to Schottky. We find that this effect be- 
comes important for inter-tube separations below 5 nm; 
the effects can be seen in the bottom curve of Fig. 2, 
where the ON state conductance saturates to a value 
much smaller than the single-tube limit because of the 
Schottky barrier. These effects of inter-tube interactions 
on contact properties can be mitigated by embedding the 
NTs in the metal contact. 
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FIG. 3. Conductance at V gs = versus NT separation, 
for different channel lengths. Solid lines are fits of the form 
A{l-e' 



To quantify the role of Coulomb interactions in the 
channel, we calculated the conductance versus gate volt- 
age for devices with different channel lengths; we plot in 
Fig. 3 the conductance at V gs — versus tube separa- 
tion, for several values of the channel length (the smallest 
value of d is large enough to avoid contact effects) . The 



solid lines in the figure are fits of the form A (l — e~ L l x ) , 
from which we extract the value of d at which the conduc- 
tance decreases by more than 10% from its d — > oo value. 
This length scale is plotted in Fig. 4 as a function of 
channel length, and delimits regions where inter-tube in- 
teractions distort the single-tube behavior. Clearly, very 
short channel devices can have very high packing densi- 
ties, while long channel devices are limited to tube sep- 
arations of 15 nm. A notable aspect of the results of 
Fig. 4 is that the value of d is independent of the chan- 
nel length for large L. At first glance, one would expect 
that a larger channel length leads to larger total charge 
on the NTs and thus larger interaction energy. However, 
screening of the Coulomb interaction by the planar gate 
leads to a different behavior, as we now discuss. 
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FIG. 4. Tube separation below which inter-tube interac- 
tions become important. Note the saturation of the length 
scale as the channel length increases. 

To derive approximate analytical expressions for the 
Coulomb interactions between NTs, we consider the sim- 
plified situation of Fig. 5a: a NT of length L is at a 
distance L g above a metallic substrate held a potential 
V g , and carries charge per unit length A. We use an im- 
age potential construction to calculate the electrostatic 
potential V at a distance d from the charged tube, giving 



\ (L + VL 2 + Ad?) -L + JL 2 + Ad 2 + 16L 2 

V = ^ + T- ln V 7 \ =4 

4lT£ (-L + VL 2 + Ad 2 ) [L + Jl 2 + Ad 2 + IQLf. 

(4) 

In the long channel limit L ^ d,L g , Eq. (4) becomes 



V = V Q + 



2tt£ 



111 



d 2 + AL 2 



(5) 



The important point here is that in this limit the poten- 
tial shift is independent of the channel length L. Taking 
a potential change of 0.1 Volts as a criterion for the im- 
portance of inter-tube effects, we obtain a separation 
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d* = 



2L n 



(6) 



below which inter-tube effects become important. Thus 
the gate oxide thickness L g sets the length scale d* . In the 
above expression, n is the effective dielectric constant for 
the device geometry. For a simple comparison with our 
simulations we take k as the dielectric constant of SiC>2, 
and use the computed value of A = 1 x 1CP 3 e/C atom 
in the ON state. This gives d* w 10 nm in good agree- 
ment with our numerical results. A key point is that the 
length scale d* depends exponentially on the gate insula- 
tor dielectric constant; thus, replacing SiC>2 with high-K 
dielectric materials (NTFETs have been fabricated with 
Hf0 2 12 , Sr0 2 13 and SrTi0 3 14 ) allows a reduction of d* 
to extremely small values, thus permitting high device 
densities. 
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FIG. 5. Sketch of situations considered for analytical calcu- 
lation of interactions between (a) semiconducting nanotubes, 
(b) semiconducting and metallic nanotubes. 

The general behavior of interactions between semicon- 
ducting NTs can be understood from further analysis of 
Eq. (5). Near threshold, the charge on the NT can be 
written as A = A (V g — Vth) where Vth is the threshold 
voltage. Substitution in Eq. (5) gives 
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(7) 



Hence, the gate voltage is rescaled and shifted by the in- 
teraction between NTs. This behavior is indicated in the 
inset of Fig. 2c, showing good collapse of the data after 
appropriate rescaling and shifting of V g . 

The results discussed above focused on interactions be- 
tween semiconducting NTs. However, in multi-NT de- 
vices, there is often a mixture of semiconducting and 
metallic nanotubes. A question therefore is how the 
presence of nearby metallic NTs affects the properties 
of the semiconducting NTs. To address this issue, we 
performed calculations for the geometry depicted in Fig. 
1(e). There, a semiconducting NT is separated from two 
metallic NTs by a distance d. Figure 6 shows the cal- 
culated conductance vs gate voltage dependence for a 
10 nm channel device. (The conductance plotted here 



is that of the semiconducting NT only. For single elec- 
trodes making contact to all the NTs, the conductance 
of the metallic tubes would have to be included to obtain 
the total device conductance.) The behavior in Fig. 6 
is similar to that observed for interactions between semi- 
conducting NTs, i.e. a degradation of the characteristics 
with decreasing separation between the NTs. However, 
because the charge on the metallic tubes is not strongly 
modulated by the gate voltage, the physics is somewhat 
different; indeed, the metallic NTs essentially consist of 
lines of constant potential. To their effect on the semi- 
conducting NTs we consider the situation of Fig. 5b, 
where a line held at a constant potential V = Vq is above 
a metallic plane held at potential V g . Solution of the 
Laplace equation for this geometry gives the potential 
on the nearby NT as 



V = V L « 



•d? 



Ll +Vg d 2 



(8) 



Thus the presence of the metallic tubes shifts and rescales 
the gate voltage seen by the semiconducting tube. For 
small separations between the NTs, the potential on the 
semiconducting NT approaches that of the metallic NT, 
which is larger than that of an isolated semiconducting 
NT. Thus, at the contact, the Fermi level is pushed above 
the valence band edge, leading to a Schottky barrier and 
saturation of the conductance to lower values. The be- 
havior expressed in Eq. (8) can be verified by looking 
at the OFF state behavior, and scaling and shifting the 
gate voltage appropriately; the inset in Fig. 6 shows the 
resulting collapse of the data. 



CD 



8 06 



cz 
o 

o 



0.4 



1 1 1 1 1 1 
. ■ No metallic tubes 1 

■ d = 10 nm 
" ■ d = 5 nm 01 


i i i i i 


■ 

1 




. — ■ — d = 3 nm 




C— 0-01 


s 

■ 






1 




-4 -2 2 4 6 8 10 

i -£ — ■— r~~ ■ i 



V (V) 

FIG. 6. Conductance of semiconducting nanotube versus 
gate voltage for a channel length of 10 nm, as a function of 
separation from neighboring metallic nanotubes. 



IV. CONCLUSION 

In summary, we have shown that Coulomb interactions 
between carbon nanotubes can have a strong influence on 
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device behavior. Below a characteristic length scale, in- 
teractions between NTs can significantly degrade device 
behavior; this can be understood in terms of a rescaling 
and shift of the applied gate voltage. For long chan- 
nel devices, the tube separation above which interactions 
between NTs is negligible becomes independent of the 
channel length, and is set by the gate oxide thickness. 
This length scale can be substantially reduced by us- 
ing high-K dielectrics, due to an exponential dependence 
on dielectric constant. While this paper focused on the 
static device properties, we expect that interactions dur- 
ing current flow (e.g. Coulomb drag) would also have 
intriguing properties. We hope that this work will stim- 
ulate controlled experiments to further explore interac- 
tions between nanotubes. 
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